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Abstract 


The lifetime of the meson is measured using the semileptonic decay 
The data sample consists of about 110 pb“^ of pp colli¬ 
sions at y/s = 1.8 TeV collected by the CDF detector at Fermilab. Four 
different decay modes are reconstructed resulting in approximately 600 
signal events. The B^ meson lifetime is determined to be t{B^) = 
(1.36 ± 0.09 loios) where the first and second uncertainties are statistical 
and systematic, respectively. The B^ meson decay length distribution is ex¬ 
amined for a lifetime difference AT/F between the two mass eigenstates of 
the B^ meson. An upper limit of AF/F < 0.83 is set at 95% confidence level. 

PACS numbers: 13.20 He, 13.25.Hw, 14.40.Nd 


Typeset using REVTgX 
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I. INTRODUCTION 


The lifetime differences between different bottom flavored hadrons probe B decay mech¬ 
anisms which are beyond the simple quark spectator model. In the case of charm mesons, 
such differences have been observed to be quite large {t{D^)/t{D^) 2.5) i. Among bot¬ 

tom hadrons, the lifetime differences are expected to be smaller due to the heavier bottom 
quark mass Some QCD inspired models based on the heavy quark expansion @] predict 
a difference between the and B^ meson lifetimes of about 5% but expect the B^ and B'^ 
lifetimes to differ by no more than 1%. Although some assumptions in Ref. have been 
questioned in Ref. |0], there is agreement that the models expect a difference between the 
B^ and B^ lifetimes of less than about 1%. These predictions are consistent with previous 
results of the B^ and B^ meson lifetimes, as well as recent B^ lifetime measurements 

In the Standard Model 0, the B^ meson exists in two UP-conjugate states, \B^) = |6s) 
and |Ps) = |5s). The two mass eigenstates of the Pg meson, B^ and B^ [H = ‘heavy’ and 
L = ‘light’), are not CP eigenstates but are mixtures of the two UP-conjugate quark states: 

l^s^) = P |i^s) - Q 14°) and \B^) = p |Ps°) + Q l4), with p^ + q^ = l. (1) 


The mass and lifetime differences between the B^ and Pg can be defined as 


Am = ruH — mi, AT = Tj 


H, 


and T = 


r^ + r, 


( 2 ) 


where mH,L and T h l denote the mass and decay width of B^ and Pg . Unlike in the case of 
the P° meson, the width difference in the B^ system is expected to be large [0 . Theoretical 


s system the 


estimates predict AT/T to be on the order of 10% to 20% ||,|]. In the Pg 
ratio Am/AT is related to the ratio of the Cabibbo-Kobayashi-Maskawa (CKM) [|^ matrix 
elements \VchVcs\l\VtsVth\-i which is quite well known, and depends only on QCD corrections 
within the Standard Model P,PT[]. Currently these QCD corrections are known to next-to- 
leading order in the l/m^ expansion §]. A measurement of AT would therefore imply a 
determination of Am and thus a way to infer the existence of Pg meson oscillations, which 
will ultimately determine the ratio of the CKM matrix elements |V)d|/|bts|- 

It is assumed that Pg mesons are produced as an equal mixture of B^ and Pg |P]. In a 
search for AT, the B^ meson decay length distribution can be to described by a function of 
the form 


J^{t) = e-^«* + e-^^* with Vl,h = T ± Ar/2, (3) 

rather than by just one exponential lifetime e“^* which is the functional form used in the 
measurement of the Pg lifetime assuming a single lifetime. 

In this paper, we present an update of the Ps lifetime measurement at CDF P| using 
the semileptonic decayQ P^ —>■ {i = e,/i), where the Pg is identified via the four 

decay modes Pg —>■ (j)TT~ , K*^K~, KgK~ , and We also examine the Pg decay length 


^Throughout the paper references to a specific charge state imply the charge-conjugate state as 
well. 
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distribution for a lifetime difference AF/F with a fit to two exponential lifetimes. The data 
sample consists of approximately 110 pb“^ of pp collisions at y^=1.8 TeV collected with 
the CDF detector during Run 1. Of this, approximately 20 pb“^ were collected during the 
1992-93 running period, while about 90 pb“^ were accumulated during the 1994-96 run 
of the Tevatron Collider. The result presented in this paper supersedes CDF’s previous 
measurement of the lifetime using semileptonic decays |^. That publication was 
based on 20 pb“^ of data and reconstructed the meson in the 07r“ decay mode only. 

The outline of this article is as follows: After a short description of the CDF detector in 
Section the selection of the Dg candidates is detailed in Section H- The determination 
of the Rs lifetime is the topic of Section We describe the search for a lifetime difference 
AF/F in Sec. 0 and offer our conclusions in Section 


II. THE CDF DETECTOR 


The Collider Detector at Fermilab (CDF) is a multi-purpose detector designed to study 
1.8 TeV pp collisions produced by the Fermilab Tevatron Collider. The detector has a coor¬ 
dinate system with the z-axis along the proton beam direction, the y-a.xis pointing vertically 
upwards, and the x-axis pointing horizontally out of the Tevatron ring. Throughout this 
article p is the azimuthal angle, 6 is the polar angle measured from the proton direction, 
and r is the radius perpendicular to the beam axis. The CDF detector is described in detail 
elsewhere [|^ . We summarize here only the detector features most relevant to this analysis. 

Three devices inside the 1.4 T solenoidal magnetic held are used for the tracking of 
charged particles: the silicon vertex detector (SVX), a set of vertex time projection chambers 
(VTX), and the central tracking chamber (CTC). The SVX [T^ consists of four layers of 


silicon microstrip detectors located at radii between 2.9 cm and 7.9 cm from the interaction 
point. It provides spatial measurements in the r-p plane with a track impact parameter 
resolution of about (13-|-40/pt) HHl; where px is the component of the track momentum p 
transverse to the 2 ;-axis (px = P ■ sin 6) given in GeV/c. The geometric acceptance of the 
SVX is about 60% as it covers only ± 25 cm from the nominal interaction point whereas the 
luminous region of the Tevatron beam has an RMS of ~ 30 cm along the beam direction. 

The VTX, which is located outside the SVX up to a radius of 22 cm, reconstructs track 
segments in the r-z plane and is used to determine the 2 ;-position of the primary interaction 
vertex with a resolution of about 0.2 cm on average. Surrounding the SVX and VTX is 
the CTC, located between radii of 30 cm and 132 cm. The CTC is a 3.2 m long cylindrical 
drift chamber that contains 84 layers of sense wires grouped into nine alternating super¬ 
layers of axial and stereo wires with a stereo angle of 3°. The outer 54 layers of the CTC 
are instrumented to record the specihc ionization dE/dx of charged particles. The CTC 
covers the pseudorapidity interval |p| less than about 1.1, where p = — ln[tan(6*/2)]. The px 
resolution of the CTC combined with the SVX is a(px)/pT = [(0.0066)^ -|- (0.0009 px)^]^'^^, 
with Px measured in GeV/c. 

Outside the solenoid are electromagnetic (CEM) and hadronic (CHA) calorimeters (|p| < 
1.1) that employ a projective tower geometry with a segmentation of Ap x Ap ~ 0.1 x 15°. 
The sampling medium is composed of scintillators layered with lead and steel absorbers. A 
layer of proportional wire chambers (CES) is located near shower maximum in the CEM and 
provides a measurement of electromagnetic shower prohles in both the <p- and 2 ;-directions. 
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The muon detection system has four of its layers of planar drift chambers (CMU) located 
beyond the central calorimeters. To reduce the probability of misidentifying penetrating 
hadrons as muon candidates in the pseudorapidity region \r]\ < 0.6, four more layers of 
chambers (CMP) are located outside the magnet return yoke. To reach these two detectors, 
particles produced at the primary interaction vertex with a polar angle of 90° must traverse 
material totaling 5.4 and 8.4 pion interaction lengths, respectively. An additional set of 
muon chambers (CMX) is located in the pseudorapidity interval 0.6 < |? 7 | < 1.0 to extend 
the polar acceptance of the muon system. 


III. THE DATA SELECTION 

In this section, we describe the data selection, which begins with the description of the 
lepton trigger data sets. This is followed by a summary of the selection requirements, which 
are applied to obtain the candidate events used for the lifetime measurements. At 

the end of this section, we briefly describe the Monte Carlo simulation of our data. 


A. The Lepton Trigger Data 


Events containing semileptonic decays are collected using inclusive electron and muon 
trigger data as well as a data set obtained from a dimuon trigger. CDF uses a three-level 
trigger system. The hrst two levels are hardware based triggers, while Level 3 is a software 
trigger based on the offline reconstruction code optimized for computational speed. 

At Level 1, inclusive electrons are selected by the presence of a single calorimeter tower 
above a threshold of 6-8 GeV depending on run conditions, while inclusive muons require 
the presence of a track in the CMU as well as the CMP. At Level 2, both of these triggers 
demand a charged track with px > 7-5 GeV/c reconstructed in the r-ip plane of the CTC 
by the central fast tracker (CFT) [|^], a hardware track processor, which uses fast timing 
information from the CTC as input. The momentum resolution of the CFT is o'(pt)/Pt ~ 
3.5% with a high efficiency. In the case of the electron trigger, this track has to be matched 
to a cluster in the electromagnetic calorimeter with transverse energy > 8.0 GeV, where 
Ft = E ■ sin^, with E being the energy of the calorimeter cluster. In the case of the muon 
trigger, this track must be matched to a reconstructed track-segment in both the GMU and 
GMP. At Level 3, a computer farm is used to fully reconstruct the data including three- 
dimensional track reconstruction in the GTG. However, the fast algorithm used for tracking 
is only efficient for particles with px > 1-4 GeV/c. In the third level of the trigger more 
stringent electron and muon selection criteria, which are similar to those described in the 
next Section |f II B| , are applied. 

The reconstruction of the —>■ decay mode is based on dimuon trigger data 

where the Level 1 trigger requires two muon candidates be observed in the muon system. 
The second level trigger requires the detection of at least one GFT track with px > 2 GeV/c 
to match in p of each muon candidate. The third level trigger requires that two reconstructed 
GTG tracks are matched with two tracks in the muon chambers and the dimuon invariant 
mass is less than 2.8 GeV/(?. During Run I about 7.5 x 10® electron trigger events and 
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about 2.5 X 10® inclusive muon trigger events were recorded by CDF, while about 1.3 x 10® 
dimuon trigger events with a dimuon invariant mass of less than 2.8 GeV/c^ were recorded. 


B. The Lepton Identification 

The identification of electron candidates reconstructed after data collection uses infor¬ 
mation from both the calorimeters and the tracking chambers. The longitudinal shower 
profile has to be consistent with an electron shower with a leakage energy from the CEM 
into the CHA of less than 4% if one track is pointing to the calorimeter tower or less than 
10% if more than one track is pointing to the calorimeter tower. The lateral shower profile 
of the CEM cluster has to be consistent with that from test beam electrons. Addition¬ 
ally, a comparison of the CES shower prohle with that of test beam electrons has to 
result in < 10. For the association of a single high px track with the calorimeter shower 
based on the position matching at the CES plane, it is required that r|A(p| <1.5 cm and 
|Azsin0| < 3 cm. In addition, we demand the Et of the electron candidate reconstructed 
offline to be greater than 6 GeV. Electron candidates from photon conversion due to detector 
material are reduced to less than 10% by looking for oppositely charged tracks which have 
a small opening angle with the electron candidate. 

The reconstruction of muon candidates is described in Ref. [|^]. We compute a 
characterizing the separation between the track segment in the muon chamber and the 
extrapolated CTC track, where the uncertainty in this variable is dominated by the 
multiple scattering in the detector material. We require < 9 in the r-ip view (CMU, 
CMP, and CMX) and < 12 in the r-z view (CMU and CMX). The transverse muon 
momentum reconstructed offline is required to be px > 6 GeV/c. For the dimuon sample 
this cut is pt > 2 GeV/c for each muon candidate. Finally, for optimal vertex resolution 
the electron and muon candidate tracks have to be reconstructed in the SVX detector. 


C. The Dg Selection 


The Ds candidates are reconstructed in the decay modes 

(i) Ds ^ 07r-, 0 ^ K+K-, 

(ii) Ds K*^K-, K*^ K+tt-, 

(hi) Ds ^ K^K-, Kl ^ TT+TT-, 

(iv) —>■ (p —>■ K^K~. 


The Dg reconstruction is explained in the next section with the example of the 07r“ 
decay channel. The other decay channels (ii)-(iv) are reconstructed in a similar way. 
We then describe only the differences in the selection of these decays (ii)-(iv) as compared 
to the Dg —> 07r“ decay mode. The kinematic selection criteria used in this analysis are 
optimized by maximizing the quantity Ns/\/Ns + Nb, where Ns is the predicted number 
of signal events based on Monte Garlo calculations (see Sec. Eni) and Nb is the observed 
number of background events estimated from the sideband regions (see Sec. [IV Ep. 
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1. The (pTT Mode 

The —»• 07r“ reconstruction starts with a search for 0 —*• K^K~ candidates. We hrst 

define a search cone around the lepton candidate with a radius AR = (Ar/)^ + {AipY of 

0.8. Any two oppositely charged tracks with px > 1-2 GeV/c within that cone are assigned 
the kaon mass and combined to form a 0 candidate. Each ej) candidate is required to have a 
mass within ±10 MeV/c^ of the world average 0 mass |^. The 0 candidate is then combined 
with another track of px > 0.8 GeV/c inside the cone which has the opposite charge of the 
lepton (we call this the ‘right-sign’ combination). This third track is assigned the pion mass. 
To ensure a good decay vertex measurement, track quality cuts requiring a minimum number 
of hits in the GTG are imposed on the tracks forming the candidate. In addition, each 
of the tracks forming the is required to be reconstructed in the SVX with hits in at 
least three out of the four silicon layers and the of fh® track fit per SVX hit has to be 
less than 6 to reject badly measured tracks. The same track selection requirements are also 
applied to the lepton candidate tracks. 

The specihc ionization information dE/dx from the GTG is used to help identify hadrons 
in the reconstruction. Because of the large Landau tail of the ionization distribution, 
the 80% truncated mean of the measured charges from the GTG sense wires is taken as the 
best estimator of the track dE/dx. The probabilities, P{i), for a track to be consistent with 
the i = e, fi, vr, K, or p hypotheses are then calculated using the measured dE/dx value and 
the predictions for the assumed particle hypotheses. We define a likelihood ratio, ^h^E/dx^ 
for a track being a kaon to be the ratio of P{K) divided by the sum of the probabilities 
of all particle hypotheses. The quantity ih'^E/dx defined correspondingly. We require the 
likelihood ratios ih/^E/dx > O-Of ^^dE/dx > 0.01. 

Since the 0 has spin 1 and both the and vr” are spin 0, the helicity angle T, which 
is the angle between the K~ and directions in the 0 rest frame, exhibits a distribution 
dN/d{cosd!) ~ cos^T. A cut |cosT| > 0.4 is therefore applied to suppress combinatorial 
background, which is found to be flat in the cos T distribution. We also apply an isolation 
cut Ex°/px(07r“) < 1.0 on the candidate, where Ef” is the sum of the transverse energy 
within a cone of radius 0.4 in rj-ip space around the lepton candidate, excluding the lepton 
energy. This cut eliminates many of the fake combinations from high track multiplicity 
jets. Finally, the mass of the D^i~^ system is required to be between 3.0 GeV/c^ and 
5.0 GeV/c^ to be consistent with coming from a decay. 

In Figure |^) the 07r“ invariant mass distribution for the ‘right-sign’ D^£~^ combinations 
is displayed. A signal of 220 ± 21 events is observed from htting a Gaussian with a 
mean of (1968 ± 1) MeV/c^ and a width of (10.5 ± 1.1) MeV/c^ plus a straight line to 
describe the combinatorial background. The mass region indicated by a dashed line has not 
been included in the fit to avoid contributions from decays where particles have been 
missed such as the 7r° in the decay > 07r“7r°. Evidence of the Gabibbo-suppressed decay 
11“ —>■ 07r“ is also present in Fig. |l|a). The shaded distribution shows 07r“ candidates from 
‘wrong-sign’ D^i~ combinations. Here, no enhancement is seen in the mass region. We 
allow multiple 07r“ combinations per event but the number of multiple entries is found to be 
less than 1% of the total number of combinations. This is also the case for the —>■ K*^K~ 
and K^K~ decay modes. 
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K*^K- Mode 


2. The D,- 


For the ^ K*^K~ decay mode with K*^ iF+vr”, we assign the and n~ masses 
to two oppositely charged tracks fonnd in the cone of A/2 < 0.8 aronnd the lepton. To 
reflect the different decay kinematics of the K*^ —>■ K^n~ decay compared to 0 —> K^K~, 
we reqnire > 1.2 GeV/c and > 0.4 GeV/c. Each K*^ candidate is reqnired 

to have a mass within ±40 MeV/c^ of the world average mass |l|. To farther rednce 
the high combinatorial backgronnd in this decay channel, we tighten the helicity cut to 
I cos4/| > 0.5 and introduce a track based isolation requirement pt(-Ds/■'■)/ pT(cone) > 0.6, 
where pT(cone) is the sum of the transverse momenta of all tracks in a cone AR < 1.0 in 
p-(f) space. The cone has the lepton direction as its axis and the primary event vertex (see 
Sec. [IV A|) as its vertex. All other selection requirements discussed in the (/)7r“ decay mode 
remain the same with the exception of the ih^E/dx described below. 

The K*^K~ invariant mass distribution for the ‘right-sign’ combinations is shown 

as dots with error bars in Figure [^d). A signal at the mass is visible. This signal 
contains also events from a D~ reflection, where the n~ is incorrectly assigned 

the kaon mass. This reflection is further discussed in Section [IV Q . To reduce the effect of 
this reflection the dE/dx requirement for that track is tightened to (-h^E/dx >0.1, while we 

still demand ^h^'E/dx ^ 0.01 for the tracks forming the K*^ candidate. 


3. The He 


K^K- Mode 


The reconstruction of the Dr 


KIK- 


decay mode begins with a search for Kg 


TT+TT” candidates by assigning the pion mass to any two oppositely charged tracks with 
Pt > 0.4 GeV/c in a cone of AR < 1.0 in rj-cj) space around the lepton. These two tracks 
are constrained to come from a common vertex and their invariant mass has to be within 
5cr of the nominal Kg mass [|l|, where a is the uncertainty on the mass measurement. 

Exploiting the long lifetime of the Kg meson, we require the Kg vertex to be significantly 
displaced from the primary event vertex, which is further described in Sec. [IV Aj . We deter¬ 
mine the transverse decay length Lxy (see Sec. [IV B|) of the Kg and require L^y > 3a, where 


a is the measured uncertainty on Lxy for each candidate event. 

The Kg candidate is combined with any kaon candidate with px > 1-2 GeV/c within 
AR < 0.8 around the lepton to find the Dg candidate. The dots with error bars in Figure ^c) 
show the KgK~ invariant mass distribution for the ‘right-sign’ DgK combinations. An 
enhancement at the Dg mass is visible. As in the K*^K~ mode, this signal contains events 
from a D~ Kg7r~ reflection, where the vr” is assigned the kaon mass (see Section Evg). 
To reduce the effect of this reflection, we again require ih^E/dx > while we demand 


ih 


dE/dx 


> 0.01 for the tracks forming the Kg. 


4- The Dg (j)p v Mode 

For the Dg —> decay mode, we start with two oppositely charged muons with 

Pt > 2 GeV/c utilizing the dimuon data set obtained with a trigger which requires the 
dimuon invariant mass to be smaller than 2.8 GeV/c^. This requirement is more than 
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90% efficient for a double semileptonic decay followed by 0/x“z/. In 

addition, two oppositely charged tracks with px > 0.8 GeV/c are assigned the kaon mass and 
combined to form a 0 candidate. There is an ambiguity in the assignment of one of the two 
muons to a found 0 — K^K~ candidate. One of the muons comes from the semileptonic 
decay while the other originates from the decay (/ieo)- ffi order to resolve this 

ambiguity, we require m.{KK^j^-) to be smaller than the world average mass |^, while 
has to be greater than m^-. To reduce combinatorial background in this decay 

channel, we use the track based isolation quantity pt(.Ds ^^)/PT(cone) and require it to be 
greater than 0.5. As required in the other decay modes, the invariant mass of the KKfifi 
system has to lie between 3.0 GeV/c^ and 5.0 GeV/c^. The number of multiple 
combinations per event is larger compared to the other three decay modes (about 10%). We 
therefore allow only one candidate per event by choosing the combination with 

the largest probability from the combined vertex £t (see Sec. rm - 

The K^K~ invariant mass distribution for the /i+ sample is shown in Figure |l|d) 
with the fit result overlaid. To obtain the number of 0 signal events, we £t a second 
order polynomial together with a Breit-Wigner line shape convoluted with a Gaussian to 
account for detector resolution. We hnd 205 ± 38 0 signal events and measure the 0 mass 
to be (1020.1 ± 0.5) MeV/c^ in agreement with the world average 0 mass |]l[. The shaded 
histogram in Fig. [^d) shows the ‘wrong-sign’ KK mass spectrum, where we consider events 
with same sign or combinations as ‘wrong-sign’. For display purposes the 

‘wrong-sign’ distribution is scaled by a factor of 0.6 to the same area as the combinatorial 
background of the ‘right-sign’ K^K~ distribution. The ‘wrong-sign’ distribution describes 
very well the shape of the combinatorial K^K~ background. No indication of a 0 signal is 
evident in the ‘wrong-sign’ distribution. 


D. The Monte Carlo Simulation 


Some quantities in this analysis like efficiencies or the iF-factor distribution further de¬ 
scribed in Sec. |1V B| are determined using a Monte Garlo (MG) calculation of b quark pro¬ 
duction and B meson decay followed by a simulation of the detector response to the hnal 
state particles. Since we extract only kinematic quantities of the B hadron decay from this 
Monte Garlo study, we do not simulate the underlying event from the pp scattering nor 
include fragmentation products, but generate only B hadrons and their decay products. 

The MG simulation begins with a model of b quark production based on a next-to-leading 
order QGD calculation |^. This calculation employs the MRSDO parton distribution func¬ 
tion to model the kinematics of the initial state partons. We generate b quarks in the 
rapidity interval |pf,| < 1.0 with a minimum px for the b quark of 8 GeV/c and 5 GeV/c to 
simulate events corresponding to the single lepton and dimuon data samples, respectively. 
These px requirements are chosen in a way to avoid any biases in the B meson kinematic 
distributions after the application of the kinematic cuts used in the analysis. The b quarks 
are fragmented into B mesons according to a model using the Peterson fragmentation func¬ 
tion with a Peterson parameter of e?, = 0.006. The bottom and charm hadrons are 
decayed into the various hnal states using branching ratios and decay kinematics governed 
by the world average masses and lifetimes of the involved particles M. 
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Events with a lepton above a momentum threshold corresponding to the appropriate 
hardware trigger are kept based on an efficiency parameterization of the CFT trigger that 
depends on the lepton px- The accepted events are passed through a simulation of the CDF 
detector that is based on parameterizations and simple models of the detector response, 
which are functions of the particle kinematics. After the simulation of the CDF detector, 
the same selection criteria applied to the data are imposed on the Monte Carlo events. 


IV. THE fiO lifetime ANALYSIS 

In this Section, we describe the measurement of the lifetime starting with the deter¬ 
mination of the primary event vertex followed by the reconstruction of the decay length. 
In order to determine the number of signal events used as a constraint in the lifetime £t, 
a reflection from D~ decays has to be considered for the —>■ K*^K~ and —> KgK~ 

decay modes. In Section [IV D| background from non-i ?2 decays is discussed, while the life¬ 
time £t is detailed in Sec. jlV E] . The B^ lifetime £t results are then presented together with 
the determination of the systematic uncertainties. 


A. The Primary Event Vertex 


The lifetime reported in this paper is based on measuring the distance between the 
primary pp event vertex and the secondary B^ decay vertex in the transverse plane. We 
hrst identify the z-position of the primary interaction vertex using the tracks reconstructed 
in the VTX detector. These tracks, when projected back to the beam axis, determine the 
longitudinal location of the primary interaction with an accuracy of about 0.2 cm along the 
beam direction. The primary vertices are distributed along the beam direction according 
to a Gaussian function with a width of ~ 30 cm. On average during Run I, the number 
of reconstructed interaction vertices in a given event follows a Poisson distribution with a 
mean of about 2.5. For the B^ lifetime measurement, we determine the 2 ;-location of the 
primary event vertex by choosing the pp interaction vertex recorded by the VTX which is 
closest to the intercept of the lepton from the semileptonic B^ decay with the beamline. 
We also require the ^-coordinates of all tracks from the decay to be within 5 cm of the 
z-coordinate of this primary vertex. 

The transverse position of the primary event vertex is determined by using the average 
beam position through the detector together with the knowledge of the longitudinal primary 
vertex position from the VTX. The average beam position is calculated offline for each data 
acquisition run. This calculation yields a transverse profile of the Tevatron beam which 
is circular with an RMS of ~ 25 pm in both the x- and p-directions. We find that the 
average beam trajectory is stable over the period that a given pp beam is stored in the 
Tevatron Collider. A detailed description of the determination of the average beamline can 
be found in Ref. For the B^ lifetime measurement, we consider only events from data 
runs with a sufficiently large number of collected events to allow a good determination of 
the run averaged beamline. In this analysis, we choose not to measure the primary vertex 
event-by-event because the presence of a second b quark decay in the event coupled with 
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the low multiplicity in semileptonic B decays can lead to a systematic bias in the lifetime 
determination. 


B. The Decay Length Reconstruction 


The tracks forming the candidate are reht with a common vertex constraint referred 
to as the tertiary vertex Vq-- The secondary vertex where the decays to a lepton 

and a (referred to as Vgo) is obtained by simultaneously intersecting the trajectory of 

the lepton track with the flight path of the candidate. Since we fully reconstruct the 
meson in the (j)7r~, K*^K~, and K^K~ decay modes, we know the flight path. In 
the —>■ channel, where we do not fully reconstruct the meson, we use the 

flight direction as a good estimate of the flight path. 

The conhdence level of the combined vertex £t is required to be greater than 1%. Fur¬ 
thermore, we require that the reconstructed decay vertex Vj^- be positively displaced 

from the primary vertex as projected along the direction of the momentum. 

The transverse decay length Lxy{B^) is dehned as the displacement X in the transverse 
plane of Vgo from the primary event vertex projected onto the momentum: 




X ■ MDSn 

Ip't(DsV+)I 


(4) 


Lxy is a signed variable which can be negative for the conhguration where the particle seems 


to decay before the point where it is produced. 


The B^ meson decay time is given by 


ct{Bl) = L 


xy 


m(i?s 

Pt{B. 


ov 


(5) 


where m(i?s) is the B^ mass [Q. Since we do not fully reconstruct the i?s meson in our 
analysis, we define the ‘pseudo-proper decay length’ 




( 6 ) 


which has a typical uncertainty of ~ 60 /xm including the contribution from the hnite size 
of the primary event vertex. In addition, we introduce a correction factor 


PT{Dsn 

Pt{BI) 


(7) 


to correct between the reconstructed and the unknown pt{B^) in the data. The 

i?s meson decay time is then given as 


Ct {B^) = Lxy 


MBs) 

PTiDsM 


X K. 


( 8 ) 


The correction between and pt{B^) is done statistically by smearing an expo¬ 

nential decay distribution with a Monte Carlo distribution of the correction factor K when 
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extracting ct{B^) from the pseudo-proper decay length in the lifetime £t as described in 
Sec. IV E . The ii'-distribution is obtained from D^i'^ combinations which originate from a 
Monte Carlo simulation (see Sec. [Ill D|) of semileptonic decays into including 

with ^ i~^X. As an example, the A'-distribution is shown for the —> 07r 

and —> 0/i“z/ decay modes in Figures and ^), respectively. The A'-distributions 

have mean values of 0.86 and 0.77 with RMS values of 0.10 and 0.12 for the —>■ 07r“ 
and modes, respectively. The A"-distribution is approximately constant as a 

function of pT{B>^i~^) for the range of pT{D^i~^) corresponding to our data. 

To ensure a precise B^ lifetime determination, we consider only B^ candidates for which 
the pseudo-proper decay length is measured with an uncertainty of less than 0.1 cm. We also 
require that the candidates have a proper decay length measured between Vgo and 
of less than 0.1 cm and that the uncertainty on this proper decay length is less than 0.1 cm. 
This requirement removes background events with very long-lived candidates, where the 
long extrapolation back to the B^ decay vertex results in a poor vertex measurement. These 
requirements have already been applied to the mass distributions shown in Fig. |^. 


C. The Determination of the Reflection from D 

The reconstructions of the decay modes into K*^K~ and KgK~ suffer from reflections 
of D~ —>■ and D~ —> A'^vr", respectively, where the tt~ is incorrectly assigned the 

kaon mass. We will discuss this reflection from D~ and the determination of the true number 
of events from the decay with the example of the —> K*^K~ mode. The effect of 
this K-7T misassignment can be seen in Figure events from a R —> D~iuX Monte Carlo 
simulation with D~ K*^n~ yield an invariant mass distribution indicated by the shape of 
the shaded area in Fig. ^c) if they are reconstructed as B^ D^iuX with K*^K~, 

misinterpreting the 7r“ as K~. A significant portion of this D~ reflection lies at the 
mass peak. 

Although we have already tightened our dE/dx likelihood ratio to better identify the K~ 
track as a kaon, CDF’s dE/dx capabilities with a tt/A' separation of about la for tracks with 
Pt greater than about 1 GeV/c are not sufficient to remove this D~ reffection. Applying a 
D~ mass veto by rejecting all K*^K~ combinations which are within a ±3a window around 
the nominal D~ mass when reconstructed as , distorts the K*^K~ mass distribution. 

It would be very difficult to estimate the remaining signal from that distribution and 
use it as input to the B^ lifetime fit. We therefore choose to measure the D~ reffection 
directly from our data and account for the D~ component in the B^ lifetime fit. We use two 
methods to determine the D~ reffection in our data. 

The first method performs a simultaneous fit to the K*^K~ and K*^n~ invariant mass 
distributions, where the mass distribution is created by switching the mass assign¬ 

ment on the K~ track to be a pion. Figure shows the K*^K~ mass distribution, while 
the corresponding mass distribution is displayed in Fig. Each distribution is 

described by a Gaussian for the corresponding D~ and signal as shown in Figures ^c) 
and plus a linear lineshape to parameterize the combinatorial background. The shape 
of the corresponding D~ or reflection as obtained from a Monte Garlo simulation is 
also included in the ht as displayed in Fig. |c) and d) as the shaded areas. The two mass 
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distributions are fit simultaneously with the number of events in the Gaussian {D~) 
signal constrained to the number of events in the corresponding D~ {D^) reflection. In 

mass values is fixed based on their nominal 
signal events and 80 ± 10 events from the 


addition, the difference between the and D~ 


mass values 


The fit returns 123 ± 25 Dc 


D reflection within the signal region dehned in Sec. [yg. The £t result is shown in 
Figure |^a)-d). We perform studies using MC pseudo-experiments to verify the validity of 
this method [^]. We find that the simultaneous fitting method returns the number of true 
Hg events in our MC studies with no bias and the error obtained from the fit to the data 
agrees with the expected uncertainty of this technique for our sample size. 

The second method for determining the amount of D~ reflection in our signal events 
exploits the difference between the D~ lifetime = (1.057 ± 0.0015) ps [|l|) and the 

lifetime = (0.467 ± 0.0017) ps |^). As described in Sec. [vg, we can determine 

the lifetime in our fit for the lifetime. We modify the fitting method used to determine 

the lifetime in the following way: We replace the exponential describing the signal 

by the sum of two exponentials, one with the lifetime and one with the lifetime (see 
Sec. [IV E| about the htting method and Sec. [V] for an example of a two-lifetime £t). We fix 
the and D~ lifetimes to their nominal values [|I| and allow the relative fractions of 

and D~ to float in the fit. With this method we obtain I 29 I 34 events and events 


‘-31 


attributed to the D reflection. We again perform studies using MC pseudo-experiments 
and verify the validity of this method to work without any bias 
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We determine the weighted average of events from both methods and obtain a 
signal of 125±20 events for the > K*^K~ decay. Both methods are also used to calculate 


the number of events and the contribution from the D reflection in the 


KIK- 


decay mode. We obtain 33 ± 8 signal events for the KgK~ mode. These numbers are 
displayed in Table | together with the numbers of signal events for the > 07 r“ and 
decay modes. As further described in Sec. jtVEj, these event numbers are used 
as a constraint in the B2 lifetime £t. 


D. The Non-Hg Background 

There are two possible sources of non-strange B meson decays which can lead to ‘right- 
sign’ D^£~^ combinations. The first process originates from the decays B^ —>■ 
and B- with the or D~^ decaying semileptonically. These decays produce 

softer and less isolated leptons than the leptons from B^ semileptonic decays. Therefore we 
expect the acceptance for this background source relative to the B^ signal to be quite small. 
We use a Monte Carlo simulation of these events and estimate their contribution Jd^d in 
the following way: 


_ _ fu + fd BR{B ^ 

^ fs BR{Bl D^ihvX) 


(9) 


We use the following branching ratios and fragmentation fractions from the Particle Data 
Group 0: BR{Bl D^^hvX) = (7.6 ± 2.4)%, BR{B D^^'>D^*^X) = (4.9 ± 1.1)%, 
fu = fd = (37.8 ± 2.2)%, and /$ = (11.2 ± 2.2)%. Crei is the ratio of efficiencies and 
acceptances for both decays obtained from a Monte Carlo simulation: 


15 












( 10 ) 


e {B ^ Di*^D^*^X) 
e{Bl^D^;kuX) ■ 


The values obtained for erei are in the order of 0.5% to 1%. The calculated fractions Jd^d for 
each decay mode are compiled in Table |. The Jd^d fraction is larger for the 
decay mode because of the on average softer B meson momentum in the dimuon data sample 
compared to the single lepton trigger events. 

The second process is a four body decay B^/B^ D^]Ci~^uX, where /C denotes any type 
of strange meson. Because of the low probability of producing ss pairs and the limited phase 
space, this process is suppressed and has not been observed experimentally [^]. Based on the 


quoted limit BR(i?°/i?+ —> D^i~^uX) < 0.9% (90% CL) and our detection efficiency 

determined from MC simulation, we expect less than 1.0% of our D^i~^ combinations to 
originate from this source. 

We also consider events from B^ X decays, with one Ds decaying semilep- 

tonically. This contribution to our B^ lifetime sample is determined from Monte Carlo 
studies in the same way as described above for the B —> D^^D^*^X background. The 
obtained fractions from these decays are small and compiled in Table |. Finally, 

backgrounds with a real meson and a fake lepton from decays such as B ^ D^^D^*^X 
with a hadron from the D^*'^ decay faking a lepton are negligible due to the low probability 
of a hadron faking a lepton. 

In summary, the contribution of all above physics backgrounds is quite small compared 
to the combinatorial background. We account for contributions from B —> D^*''X and 

Bg X decays in our lifetime £t described next, and treat the contribution 

of B^/B^ D^lCt^vX decays as a source of systematic uncertainty in the i?s lifetime 
measurement. 


E. The Description of the Lifetime Fit 

As input to the B^ lifetime fit, we define a signal sample using a mass window 
from 1.944 GeV/c^ to 1.994 GeV/c^ for the —»• K*^K~, and KgK~ decay modes 

and a 0 signal window from 1.0094 GeV/c^ to 1.0294 GeV/c^ for the —>■ decay 

channel. The numbers of events in the signal samples can be found for the four decay modes 
in Table |. To model the pseudo-proper decay length distribution of the combinatorial 
background events contained in the signal sample, we define a background sample which 
consists of ‘right-sign’ events from the sidebands (1.884 GeV/c^ - 1.934 GeV/c^ and 
2.004 GeV/c^ - 2.054 GeV/c^) and ‘wrong-sign’ events from the interval 1.884 GeV/c^ 
to 2.054 GeV/c^. For the 0/i“z/ decay mode the 0 sidebands are dehned from 

0.9844 GeV/c^ to 1.0044 GeV/c^ and from 1.0344 GeV/c^ to 1.0544 GeV/c^, while the 
‘wrong-sign’ combinations are taken from the region 0.9844 GeV/c^ to 1.0544 GeV/c^. We 
assume the combinatorial background to originate from random track combinations and 
therefore use the sidebands to model the background in the signal sample. This assumption is 
supported by the mass distribution of the ‘wrong-sign’ combinations where no enhancement 
is visible at the mass. By adding the ‘wrong-sign’ combinations to the ‘right-sign’ 
sideband events, we better constrain the shape of the combinatorial background events 
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in the signal samples for decay channels with low combinatorial background like the 
> 07r“ mode. 

The pseudo-proper decay length distribution obtained from the signal sample is £t using 
an unbinned maximum log-likelihood method. Both the lifetime, denoted as cr be¬ 
low, and the background shape are determined in a simultaneous £t using the signal and 
background samples. Thus the likelihood function £ is a combination of two parts: 

Ns Nb 

£ = n + (1 - -11^4. (11) 

* j 


where Ns and Nb are the number of events in the signal and background samples, fsig is the 
ratio of signal events obtained from the mass distributions (see Table |) to the total 
number of events in the signal sample. To constrain /gig we factor in an additional term 
to the likelihood function C above to constrain the number of signal events obtained 
from the invariant mass distributions within their uncertainty. 

The signal probability function Bsig consists of a normalized decay exponential function 
convoluted with a Gaussian resolution function Q and is smeared with a normalized K- 
distribution H{K): 


= I dKn{K) 


K K T 

— exp{-} (g) ^(A*| X, sa*) 

Lcr CT 


( 12 ) 


Here, A* is the measured pseudo-proper decay length with uncertainty ex* and x is the true 
pseudo-proper decay length. Because of systematic uncertainties in the overall scale of the 
decay length uncertainties, which we estimate on an event-by-event basis, we introduce a 
scale factor, s, which is a free parameter in the lifetime £t. We subsequently vary s in 
the hts to determine the sensitivity of our measurement to this uncertainty. The integration 
over the momentum ratio K is approximated by a hnite sum 


J dKn{K) ^ AKHiKi), (13) 

i 

where the sum is taken over bin i of a histogrammed distribution T-C{Ki) with bin width AK 
as shown e.g. in Figure 

The background probability function B^g is parameterized by a Gaussian centered at zero, 
a negative exponential tail, and a positive decay exponential to characterize the contribution 
of heavy flavor decays in the background sample: 

B^^{x) = {l-U-f.)g{N\x,sa^) + 

+ ^exp{--^} 0 g{X"\ X, S(T®) -F 

A+ A+ 

-h ^exp{--^} 0 ^(A*| X, S(J*). (14) 

A_ A_ 


Here, f± are the fractions of positive and negative lifetime backgrounds and A± are the 
effective lifetimes of those backgrounds. We verify the parameters f± and A± agree with 
the ‘right-sign’ sideband events and the ‘wrong-sign’ combinations separately, allowing us 
to combine both samples resulting in the background samples described above. 
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The events originating from the D reflection in the K*^K and 


decays (see Sec. |IV C|) are also accounted for in the likelihood function by a term 


KIK- 


dKn{K) 


fo- 


K 


ct{B) 


exp{- 


Kx 


ct{B) 


} ® ^(A*| X, scr* 


(15) 


where fjy- refers to the fraction of the D reflection in the sample and ct{B) is taken 
to be the world average B^ lifetime [|l|. 


F. The Fit Results 

We hrst determine the B^ lifetime for each of the four decay channels individually. 
The parameters allowed to float in the fit are the B^ lifetime, /sig, A±, f±, and the overall scale 
factor s. The htted values for ct{B^) and their statistical uncertainties are shown in Table |, 
and are in good statistical agreement. The pseudo-proper decay length distribution of the 
signal sample with the result of the £t superimposed is shown in Figure^) for the —»• 07 r“ 

decay mode. The dashed line represents the B^ signal contribution, while the shaded curve 
shows the sum of the background probability function over the events in the signal sample. 
The same distribution of the background sample is displayed in Figure ij) with the result 
of the fit superimposed. Figures & and 1 ^ show the corresponding distributions for the 
Ds ^ K*^K-, Ds K^K-, and v decay modes, respectively. The combined 

B^ lifetime from all four decay modes is determined from a simultaneous £t to be 
ct(B^) = (408 I 27 ) or t(B^) = (1.36 ± 0.09) ps, where the errors shown are statistical 
only. 

As a consistency check, we use the —>■ ((> 7 r“ decay mode to also £t the B^ lifetime 
from the proper decay length measured from the secondary vertex Vb^ to the tertiary ver¬ 
tex Vd^- Since the B^ decay is fully reconstructed, its relativistic boost is known and a 
convolution with a px-correction factor distribution in the fit does not apply. The result is 
ct{BI) = (136 IJ 5 ) pm (statistical error only), which is consistent with the world average 
B^ lifetime [|[]. Figures ||a) and ||b) show the B^ proper decay length distributions for the 
signal and background samples, respectively, with the results of the fit superimposed. 


G. The Systematic Uncertainties 


Table |n| lists all sources of systematic uncertainty considered in this analysis. The major 
contribution originates from the treatment of the background. In particular, the following 
evaluations of systematic errors are performed yielding the uncertainties reported in Tab. ||: 


Background treatment: The combinatorial background in the signal sample is param¬ 
eterized by the positive and negative lifetimes A_|_ and A_ as well as their respective 
fractions /+ and /_ as described in Eq. ([T^). To evaluate the dependence of the 
Bg lifetime on the background parameterization, we vary A± and /± within Tier of 
their values returned from the B^ lifetime £t. Since there is a correlation between the 
background parameters \± and f± and the B^ lifetime in the simultaneous £t to the sig¬ 
nal and background samples, part of this systematic uncertainty is already accounted 
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for in the statistical error on However, this correlation is small. We therefore 

adopt this method as a conservative way to evaluate the systematic uncertainty from 
background treatment. 


• Non-Hg backgrounds: These backgrounds (see Sec. lYD enter our £t as hxed frac¬ 

tions. We vary these fractions within ±50% of their reported values (see Tab. ||) and 
repeat the £t. We also evaluate the background from ICi'^uX decays 

and consider its fraction to be 2.0%, twice the quoted limit [|I], pO|] . 

• Decay length requirement: There are two requirements that can bias the B^ lifetime 
result. These are the cut on |cr(Ds )| <0.1 cm and the requirement that the recon¬ 
structed decay vertex be positively displaced from the primary event vertex. 
To study the effect of these cuts, we use high statistics Monte Carlo samples. We 
first £t the lifetime with all the selection requirements, and then remove each cut 
individually noting the shift in the B^ lifetime. 


Momentum estimate: The B^ lifetime result is sensitive to the distribution of the 
correction factor 7i(/C), which can be affected by the lepton cut and the decay 
kinematics. For the standard £t we require Pt(/^) > 2.0 GeV/c for the —>• 0/r“z/ 
decay mode and Pt(^) > 6.0 GeV/c for the other decay channels. To test the effect 
of the lepton pt dependence, we generate new iF-distributions for lower and higher 
lepton Pt cut values. We also compare the effect on the kinematics of semileptonic 
decays using a pure V-A decay versus semileptonic decays using the ISGW form 
pT| . In addition, an alternative pt spectral shape of b quark production is con- 


factor 

sidered based on a comparison of the lepton pt shape in the data and in Monte Garlo 
events. Finally, the iF-distribution is somewhat dependent on the electron identifica¬ 
tion. We study a possible incompleteness in the treatment of the electron selection 
with our Monte Garlo simulation and assign a systematic error of ±3 pm. 


• Decay length resolution: Our uncertainty on the estimate of the decay length resolution 
is expressed in the scale factor s, which is htted to 1.29 ± 0.03. We £x the scale factor 
at 1.0, and again at 1.38, the latter corresponding to a ±3cr upward shift from the 
htted value, and repeat the B^ lifetime £t. 

• D~ rehection: The rehection from D~ in the K*^K~ and decay 

modes changes the number of signal events in these two channels. We study the 
inhuence of the D~ rehection by varying the number of signal events within their 
error as determined in Sec. 


Detector alignment: We also account for a possible residual misalignment of the SVX 
and assign an error of ±2 pm as further detailed in Ref. |]T^ . 


The systematic uncertainties noted above have been combined in quadrature. Quoting 
the statistical and systematic uncertainties separately, we measure the B^ lifetime using 
semileptonic B^ decays to be 


r(B”) = (1.36 ± 0.09 PS. 


(16) 
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where the hrst error is statistical and the second systematic. This resnlt is currently the 
world’s best measurement of the lifetime from a single experiment. In comparison, 
the world average lifetime is (1.57 ± 0.08) ps [^. This measurement supersedes CDF’s 
previously published B^ lifetime result of t{B^) = (1.42 ±0.11) ps using a data sample 

corresponding to an integrated luminosity of 20 pb“^ p. 

Using the CDF average B^ lifetime t{B^) = (1.513 ± 0.053) ps we determine the 
B^/B^ lifetime ratio to be 0.899 ± 0.072 taking correlated systematic uncertainties into 
account. However, ignoring the correlated systematic uncertainties increases the error on 
the B^/B^ lifetime ratio only to ±0.077, since t{B^)/t{B^) is dominated by the statistical 
error on the B9 lifetime measurement. 


V. THE DETERMINATION OF AF/F 


We examine the B^ meson pseudo-proper decay length distribution from D^i~^ correla¬ 
tions for a lifetime difference AF/F between the two mass eigenstates of the B^ meson, B^ 
and i?s • In the case of a lifetime difference in the B^ system, the decay length distribution 
for events from the semileptonic decay B^ is expected to be governed by the 

sum of two exponentials. We expand the likelihood £t to describe the B^ pseudo-proper 
decay length distribution to a function of the form 

AF 1 AF 

X{t) = e-"-* ± e-"-* with F^,^ = r± — = F-(1±-—), (17) 


rather than htting for just one exponential lifetime e”'"^ The parameter AF/F is the pa¬ 
rameter we £t for. Since AF/F is symmetric about zero, it is required to be positive. In 
the case of a lifetime difference AF ^ 0, the total decay width F = l/2-(Fi^±Fi) and the 
mean B^ lifetime rm(i?s) obtained from a fit assuming a single B^ lifetime, are no longer 
reciprocal to each other but follow the relation 


TmiBl) 


1 1 ±(#)^ 

F l-(#)2- 


(18) 


We incorporate the relation in Eq. (|^) into our likelihood htting function. We follow the 
suggestion given in Ref. 0 and £x the mean B^ lifetime to the world average B^ lifetime 
since both lifetimes are expected to agree within 1% |^,^]. This theoretical assumption can 
be verihed by the current world average t{B^) = 1.55 ± 0.05 ps and t{B^) = 1.57 ± 0.08 ps. 
The CDF average B^/B^ lifetime ratio derived above also supports this assumption. 

The £t returns AF/F = 0.34 where the given error is statistical only. This indicates 
that with the current statistics of our D^i~^ sample we are not sensitive to a B^ lifetime 
difference. Based on this fit result, we integrate the normalized likelihood as a function of 
AF/F and hnd the 95% conhdence level (CL) limit at 


AF 


< 0.83 


(95% CL). 


( 19 ) 


This is the first experimental result for the lifetime difference in the B^ system. 
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Using a value of AF/Am = (5.6 ± 2.6) ■ 10“^ from Ref. 0 and setting rm(Ss) to the 
world average lifetime |l[], an upper limit on the Rg mixing frequency of Ams < 96 ps“^ 
(95% CL) can be determined within the Standard Model. Including the dependence on 
AF/Am and in our limit, we obtain 


Arris < 96 ps ^ x 


/5.6 ■ 10-3\ 
V AF/Am ; 


X 


/1.55 ps\ 


(95% CL). 


( 20 ) 


VI. CONCLUSION 


We have presented a measurement of the B^ meson lifetime using semileptonic B^ decays, 
where the meson is reconstructed through the four decay modes B^ —> 

Bs K^K-, and B^ We obtain 

t(B“) = (1.36 ± 0.09 tS) ps. (21) 


where the hrst error is statistical and the second systematic. This is currently the world’s best 
measurement of the B^ lifetime from a single experiment. This result agrees with an earlier 
CDF measurement |^, which is superseded by the present measurement. We determine the 
B^/B^ lifetime ratio to be 0.899 ± 0.072 using the CDF average B^ lifetime . 

In addition, we have examined the B^ meson pseudo-proper decay length distribution 
for a lifetime difference AF/F between the two mass eigenstates of the B^ meson, B^ 
and Rg. Using all four B^ decay modes, an upper limit of AF/F < 0.83 is set at 95% CL, 
corresponding to the Standard Model limit 


Ams < 96 ps 


-1 


X 


/5.6- 10-3' 
i AF/ Am 


X 


1.55 




MBDJ 


(95% CL). 


( 22 ) 


With considerably increased statistics in the next run of the Tevatron Collider, our sensitivity 
on the lifetime difference AF/F will be signihcantly improved 
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TABLES 


Dg Decay Mode 

NiDs) 

JdsD 

/TsTs 

Aevt 

cr(B^) 

4>tt~ 

220 ± 21 

2.6% 

0.8% 

350 

41811^ 

K*^K- 

125 ± 20 

2.5% 

0.8% 

820 

411 til Mm 

KiK- 

33 ±8 

1.8% 

0.6% 

146 

397 till Mm 

(f)^ V 

205 ± 38 

5.7% 

1.7% 

635 

399 urn 


TABLE I. Summary of results for the four decays: 

- the number N(D^) of signal events as input to the lifetime fit; 

- the expected fraction fn^D oi B ^ decays; 

- the expected fraction Jd^d^ of decays; 

- the number A"evt of events in the signal samples; 

- the htted B^ lifetimes ct(B^), where the errors shown are statistical only. 


Error Source 

Acr(BO) 

Background treatment 

±11 filQ 

Non-Bg backgrounds 

±6 pm 

Decay length requirement 

Momentum estimate 

tl pm 

Lepton pt dependence 

Is Mm 

B decay model 

1? pm 

b quark px spectrum 

±5 pm 

Electron selection 

±3 pm 

Decay length resolution 

+7 

-2 

D~ reflection 

±1 pm 

Detector alignment 

±2 pm 

Total 

-il Mm 


TABLE II. Compilation of systematic uncertainties in the measurement of the B^ lifetime 
combining all four decay modes. 


24 












Candidates / 10 MeV/c Candidates / 10 MeV/c 


FIGURES 





FIG. 1. Invariant mass distributions of (a) , (b) —> K*^K~, (c) —> KgK~, 

and (d) cj) —> K^K~ from (j)fj,~u. The dots with error bars are for ‘right-sign’ D^£~^ 

combinations while the shaded histograms show the corresponding ‘wrong-sign’ distributions. In 
(a) evidence of the decay D~ (pTr~ is present. The results of the fits described in the text are 
also superimposed. The mass regions indicated by a dashed line have not been included in the fits. 
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Pt{Ds£+)/pt{BI) 


Pt{Ds£+)/pt{BI) 


FIG. 2. Normalized iF-factor distributions for Monte Carlo 

decays with (a) D7 (/)7r“ and (b) D7 (l)p~v. 





Entries Candidates /10 MeV/c 



FIG. 3. (a) Mass distributions for candidates in the —> K*^K~ decay mode, (b) Mass distri¬ 
bution if these candidates are assumed to be D~ —> . (c) Distribution of the K*^K~ 

signal and the reflection from D~ —> (shaded area) as obtained from Monte Carlo simu¬ 

lations. Normalizations are determined from the simultaneous fit described in the text, (d) Mass 
distribution of the corresponding D~ —> K*^7r~ signal and the reflection from —> K*^K~ 

(shaded area). 
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Candidates / 0.01 cm Candidates / 0.01 cm 



Pseudo-Proper Decay Length [cm] 



Pseudo-Proper Decay Length [cm] 


FIG. 4. (a) i?s pseudo-proper decay length distribution for the > 4>tt~ signal sample with 

the result of the fit superimposed. The dashed line is the signal contribution, while the shaded 
curve represents the contribution from the combinatorial background, (b) Pseudo-proper decay 
length distribution for the background sample with the fit result superimposed. 
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Candidates / 0.01 cm Candidates / 0.01 cm 



Pseudo-Proper Decay Length [cm] 


10 ^ 

10 ^ 

10 
1 

-0.1 0 0.1 0.2 0.3 

Pseudo-Proper Decay Length [cm] 



FIG. 5. (a) i?s pseudo-proper decay length distribution for the —> K*^K~ signal sample 

with the result of the fit superimposed. The dashed line is the signal contribution, while the 
shaded curve represents the contribution from the combinatorial background, (b) Pseudo-proper 
decay length distribution for the background sample with the fit result superimposed. 
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Candidates / 0.01 cm Candidates / 0.01 cm 



Pseudo-Proper Decay Length [cm] 



Pseudo-Proper Decay Length [cm] 


FIG. 6. (a) i?s pseudo-proper decay length distribution for the —> KgK~ signal sample 

with the result of the fit superimposed. The dashed line is the signal contribution, while the 
shaded curve represents the contribution from the combinatorial background, (b) Pseudo-proper 
decay length distribution for the background sample with the fit result superimposed. 
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Candidates / 0.01 cm Candidates / 0.01 cm 



Pseudo-Proper Decay Length [cm] 



Pseudo-Proper Decay Length [cm] 


FIG. 7. (a) i?s pseudo-proper decay length distribution for the signal sample 

with the result of the fit superimposed. The dashed line is the signal contribution, while the 
shaded curve represents the contribution from the combinatorial background, (b) Pseudo-proper 
decay length distribution for the background sample with the fit result superimposed. 
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Candidates / 0.005 cm Candidates / 0.005 cm 



Proper Decay Length [cm] 



Proper Decay Length [cm] 


FIG. 8. (a) proper decay length distribution ct{D^) for the —> 4>tt~ signal sample 

with the result of the fit superimposed. The dashed line is the signal contribution, while the 
shaded curve represents the contribution from the combinatorial background, (b) Proper decay 
length distribution ct{D^) for the background sample with the fit result superimposed. 
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